Submitted to Phys. Rev. B.(Dec 2008) 
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A metamaterials-based approach to making a wide-angle absorber of infrared radiation is de- 
scribed. The technique is based on an anisotropic Perfectly Impedance Matched Negative Index 
Material (PIMNIM) . It is shown analytically that a sub- wavelength in all three dimensions PIMNIM 
enables absorption of close to 100% for incidence angles up to 45 deg to the normal. A specific im- 
plementation of such frequency-tunable PIMNIM based on plasmonic metamaterials is presented. 
Applications to infrared imaging and coherent thermal sources are described. 



PACS numbers: 

I. INTRODUCTION 



The emergence of a new field of electromagnetic meta- 
materials was brought about by the demand for materi- 
als with exotic properties unattainable in nature. One 
such property is a negative refractive index that requires 
both the effective dielectric permittivity e e g and mag- 
netic permeability /i e ff to be negative^. Applications of 
negative index metamaterials (NIMs) include "perfect" 
lenses, sub- wavelength transmission lines and resonators, 
miniature antennas, among others— i 4 ^. 

The results reported in this paper pertain to two other 
recently emerged applications of metamaterials. The first 
one is wavelength-selective infrared^ and terahertz^ de- 
tection, important for thermal imaging, night vision, and 
nondestructive detection. Wide-angle power absorption 
efficiency is desirable for miniaturizing photodetectors 
or microbolomctcrs down to the wavelength size. The 
second application is the development of coherent ther- 
mal emitters ^ 10 ' 11 ' 12 for spectroscopic and thermophoto- 
voltaic (TPV) 13 i 14 applications. By virtue of the Kirch- 
hoff's Law, emissivity of a thermal emitter approaches 
the blackbody limit only if the absorptivity approaches 
unity. Moreover, wavelength-selective radiators can dra- 
matically improve the efficiency of current generation in a 
TPV cel l 13 ' 14 if their emission spectrum is matched to the 
bandgap of the TPV converter. For example, a typical 
TPV converter, GaSb, has the bandgap of Eq = 0.7 eV 
that would be ideally suited to a wavelength-selective ra- 
diator operating in near-infrared (around A = 1.7/xm). 
For spectroscopic applications, arrays of wavelength- 
sized narrow-band coherent thermal emitters based on a 
metamaterial unit cell can be designed. The emitted ra- 
diation of each detector is focusable to wavelength-sized 
spots using the far-field by large numerical aperture op- 
tics. With these applications in mind, we describe in 
this paper a ultra-thin metamaterials-based wide-angle 
absorber of near-infrared radiation. The design of this 
perfect absorber is inspired by a Perfectly Impedance 
Matched Negative Index Metamaterial (PIMNIM) con- 
structed from plasmonic wires. 



II. THEORETICAL BACKGROUND 

As a background, we consider a simple problem of ra- 
diation absorption by a semi-infinite slab of a lossy meta- 
material with engineered dielectric permittivity and mag- 
netic permeability tensors e and jl. Radiation is assumed 
to be incident in the x — z plane at an angle 9 with respect 
to the vacuum-material interface normal e x . We further 
assume s— polarization of the incident wave as shown in 
the inset of Fig. [I] i.e., the only non- vanishing field com- 
ponents are E y , H x , H z . Assuming that both e and p, are 
diagonal tensors, their only relevant components are e yy , 
fi xx , and /i zz . For our semi- infinite slab (assuming that 
the metamaterial's thickness L x is sufficient to absorb 
all transmitted radiation), absorptivity A is limited only 
by reflection . A straightforward calculation yields the 
reflection r and absorption A = 1 — |r| 2 coefficients: 



cos 9 - yj e yy /^i zz - sin 2 9/ (/x zz /i xx ) 

A = 1 - , (1) 

cos 9 + yj e yy /fi zz - sin 2 9/ (ji zz fJ, xx ) 

where all material parameters are, in general, 
wavelength-dependent. 

Studying Eq. (P), we note that total absorption 
(A(Xq) = 1) at a specific wavelength Ao is achievable 
at normal incidence (9 = 0) if the absorber material's 
impedance r\ = *J ^ zz /e yy is matched to that of vacuum, 
i.e., when /i zz (Ao) = e yy {\o). Using impedance match- 
ing to achieve total absorption is a very well known ap- 
proach in microwave enginccrin g 15 i 16 . More surprising 
is another prediction of Eq. (JTJ) : that nearly total ab- 
sorption can be achieved over a very broad range of an- 
gles. Assuming that ^l xx {\q) = 1 and \\i zz \ = \e yy \ ^> 1, 
we find that A(A ) » 1 - tan 4 0/2. For a 90 deg full- 
angle (9 ma x = 7r/4) we find that 0.97 < A < 1 for 
< 9 < 9 rnax . Even more remarkably, a similarly high 
broad-angle absorption is predicted even for much mod- 
est values of /i zz and e yy . For example, if the absorb- 
ing medium is a Negative Index Material (NIM}^! 2 . with 
e yy = pizz = — 1 + i, we find that 0.94 < A < 1. The 
angular dependence A(9) for such perfectly impedance 
matched negative index material (PIMNIM) is shown in 



2 



Fig. [T]for the semi-infinite metamaterial slab, L x = oo, 
as well as for the more general case of a finite slab, with 
L x ~ Ao/7. The two curves are barely distinguishable 
because of the high metamaterial loss. 
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FIG. 1: (Color online) Angular dependence of the absorption 
coefficient for the idealized structure with [i zz = e yy = — 1-H, 
Hxx = 1, and L x = oo (dashed line), L x — 200 nm (dashed- 
dotted line), barely distinguishable from the L x — oo line. 
Solid line: same for the specific implementation of the meta- 
material absorber shown in Fig. [2] A schematic of the ideal- 
ized wide-angle metamaterial absorber is given in the inset, 
S-polarized incident wave is assumed. 



The main implication of these results is that a sub- 
wavelength slab of almost any impedance-matched meta- 
material acts as a wide-angle absorber. The challenge 
now is to design such a metamaterial, especially in the 
technologically important infrared part of the spectrum. 
With this goal in mind, the rest of this paper is orga- 
nized as follows. In section UTT1 we present one such pos- 
sible metamaterial design based on two layers of alter- 
nating long and short plasmonic nanoantennas (Fig. [5]), 
and demonstrate that the unit cell of this PIMNIM can 
be made highly sub-wavelength in near-infrared. The 
sub- wavelength requirement is critical for achieving wide- 
angle absorption: recent numerical simulations^ have 
found that A(9) rapidly drops with 9 when the unit 
cell is too large. Then, in section IIV1 we demonstrate 
that a single layer of PIMNIM enables optical absorp- 
tion 0.7 < A(9) < 0.9 for a 90deg full-angle scan, and 
that its small-angle absorption coefficient exactly agrees 
with the prediction of Eq. . Finally, in section [V] we 
show that extremely high optical intensity enhancements 
(at least xlO 3 ) can be produced at the absorption peak. 



III. PIMNIM DESIGN 

The geometry of the unit cell of the proposed PIMNIM 
is shown in Fig. [5] Each unit cell consists of two par- 
allel layers separated by the distance h x . Each layer 



consists of a cut-wire of width w z and length w y sur- 
rounded by the continuous in the y— direction wires of 
width d z . To enable planar fabrication, the heights t x 
of the cut and continuous wires are assumed to be the 
same. Metal wires are assumed to be embedded inside a 
dielectric with = 2.25. An s-polarized electromagnetic 
wave E \\ e y excites the electric and magnetic responses 
of the unit cell. The magnetic response is caused by 
counter-propagating currents flowing through the adja- 
cent cut- wires. Both cut and continuous wires contribute 
to the electric response. A similar structure has been an- 
alyzed in the microwave part of the spectrum^, where it 
is not sub-wavelength. Plasmonic effects, i.e. taking into 
account that metals have a finite dielectric permittivity 
e m = e' m + ie" t with e' m < 0, are necessary to miniaturize 
the unit cell. 
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FIG. 2: (Color online) Schematic of the PIMNIM structure. 
Unit cell for electromagnetic and electrostatic simulations is 
inside the dashed rectangles. 



Electromagnetic resonances in effective permittivity 
e e ff and magnetic permeability /z e ff of plasmonic com- 
posites are unambiguously related 18 i 19 t 20 to the elec- 
trostatic surface plasmon resonances of the appropriate 
symmetry (electric dipole and magnetic dipole, corre- 
spondingly). Therefore, the first step in investigating 
the suitability of a particular unit cell geometry for a 
plasmonic NIM is to identify the frequencies of its elec- 
trostatic resonances. Specifically, magnetic activity (in- 
cluding negative index behavior) has been shown to ex- 
ist^ only for A > A cs , where A cs = 2irc/u! cs is the vac- 
uum wavelength corresponding to the frequency of the 
electrostatic resonance responsible for the magnetic ac- 
tivity. Such magnetically-active (MA) resonance of the 
PIMNIM structure has been calculated by solving the 

Poisson's equation V- ^eV^j = for the electrostatic po- 
tential 0, where e is a piecewise constant function equal 
to e m inside the metal wires and outside. Poisson's 
equation can be solved as a generalized eigenvalue equa- 
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tio n 20 i 21 i 22 to yield the values of e m corresponding to the 
electrostatic resonances. 

Potential distribution of the lowest-frequency MA res- 
onance is shown in Fig. [3] Its magnetic nature can be 
deduced from the electric field loops formed between the 
two cut-wires. The above resonance occurs when the di- 
electric contrast between metal wires and the dielectric 
matrix is e m /ed — —26.8. Because the resonant dielec- 
tric contrast is determined from electrostatic calculations 
that contains no spatial scale, this value is determined by 
the geometric shape of the unit cell and its inclusions but 
not by the overall scale set by either of the periods L XjVtZ . 
Assuming gold as the plasmonic component and silica 
with td — 2.25 as an embedding dielectric, the electro- 
static resonance occurs at A cs = 1.2 /im, thereby setting 
the lower limit on the wavelength at which strong mag- 
netic response is expected. The next, electrically active 
(EA) resonance corresponds to e m /ed = —19, is consid- 
erably blue-shifted from the MA resonance. Therefore, 
the PIMNIM structure depicted in Fig. [3] is promising 
as a NIM that has widely separated electric and mag- 
netic resonances. Thus, its electromagnetic character- 
istics strongly resemble those of the original microwave 
NIMs 2 . The important difference is that the proposed 
structure operates at the near-infrared frequencies and 
exploits plasmonic resonances to achieve sub- Ad cell size, 
where Ad is the wavelength inside the substrate. 




FIG. 3: (Color online) Electrostatic magnetically-active plas- 
mon resonance of the PIMNIM structure shown in Fig. [2]with 
the scalable geometric proportions corresponding to (in nm): 
L y = 320, w y = 200, t x = 80, w z = 50, d z = 80, h x = 20, 
v x = 10, h z = 20. The resonance occurs at e m /ed = —26.8. 
Arrows: electric field; color: electrostatic potential. 

It is instructive to compare the resonant frequency of 
the MA resonance of the PIMNIM structure with that 
of the traditional Double Fishnet (DF) 23,24 structure. 
DF can be obtained from the PIMNIM by extending 
the cut-wires in z-direction until they merge with con- 
tinuous wires running along the y-axis, thereby forming 
the second, orthogonal set of continuous wires in the z- 
direction. For example, the DF structure thus obtained 



from the PIMNIM structure shown in Fig. [3J is periodic 
in the y — z plane with the periods of L z — 170 nm and 
L y = 320 nm containing two sets of intersecting metal 
strips in y and z direction, with the corresponding widths 
of d z = 80 nm and w y = 200 nm. The MA plasmonic res- 
onance is found at e m /^d = —8.3 which corresponds to 
Acs = 0.73 fim for gold. Of course, the physical sizes 
can be scaled from the above dimensions by an arbitrary 
factor because there is no physical scale in electrostat- 
ics. This simulation illustrates the challenge of making a 
strongly sub- Ad DF structure. For the unit cell to be sub- 
Ad, the operational wavelength must be reasonably close 
to A cs . Therefore, making a unit cell with the largest 
dimension of A/4 requires L y < 200 nm. Such small DF 
structures have never been fabricated to date, which ex- 
plains why sub-Ad DF-based NIM have never been pro- 
duced. 



IV. SIMULATION RESULTS 

Full electromagnetic finite elements frequency do- 
main (FEFD) simulations of light transmission/reflection 
through the PIMNIM structure shown in Fig. [2] (ap- 
proximate physical parameters indicated in Fig. [3] and 
further optimized to achieve perfect impedance match- 
ing) were carried out using the COMSOL 25 commer- 
cial package. Drude model for e m (cj) of gold, with 
lu p = 1.367 x 10 16 rad/sec and 7 = 4 x 10 13 rad/sec 
(taken from 2 ^), and ed = 2.25 were used in the simu- 
lations. Reflection minimization was carried out using 
standard non-linear multivariable optimization to fine- 
tune the PIMNIM parameters and adjust the perfect- 
matching wavelength Ao- The normal incidence (8 = 0) 
results are shown in Fig. 0Ja): zero reflection is achieved 
(by design) at A = 1.5 /zm, with the single PIMNIM 
layer absorption coefficient A w 0.9 (or A w 0.99 for two 
layers) 

Effective parameters e yy = e c ff(A) and [i zz = ^ g(X) 
of the PIMNIM structure were extracted from the com- 
plex transmission and reflection coefficients through a 
single PIMNIM layer using the standard procedure^!. 
Because the negative index band is the lowest (in fre- 
quency) transmission band of the PIMNIM, there is no 
ambiguity in determining e c g and /^ e ff ■ At the impedance 
matching point Ao = 1.5 /im it is found that e c g = 
fi e S = —0.85 + lAi. Although the so-called figure-of- 
merit (FOM=Re n/Im n) is less than unity, achieving 
high FOM is not necessary (or even desirable) for accom- 
plishing total wide-angle absorption. More important is 
ensuring that the unit cell of the PIMNIM is very sub- 
wavelength: L x _ y _ z <C Xd- This is indeed accomplished 
by the present PIMNIM design: n d L y /\ « 0.3. 

The peak absorption, A max (6 l ) = 1 — H 2 (0), of the 
PIMNIM is plotted in Fig. [TJ For small 8 there is an ex- 
cellent agreement between A max and Eq. Jl} indicating 
that the PIMNIM is accurately described as an effective 
medium. Even for 8 = ±45 deg absorptivity remains at 
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FIG. 4: (Color online) (a) Reflectance, transmittance and ab- 
sorbance at normal incidence of a single PIMNIM layer, (b) 
Extracted effective dielectric permittivity e e ff and magnetic 
permeability fi e g of the PIMNIM structure. Reflectance van- 
ishes at A = 1.5 /im because e a s = fi c f{. Parameters (in nm): 
L y = 326, w y = 208, t x = 80, w z = 46, d z = 78, h x = 21, 
h, = 20. 



about 70%. This implies that, in a narrow spectral inter- 
val around Ao , emissivity of a heated PIMNIM is close to 
that of a blackbody over a broad range of angles, which 
would enable to manufacture wavelength-size coherent 
thermal emitters. 



V. FIELD ENHANCEMENT 

Physically, the absorption peak corresponds to strong 
electric field enhancement inside the PIMNIM struc- 
ture. For example, for the PIMNIM described by Fig. |H 
it is found that the maximum intensity enhancement 
near the corners and edges of the cut-wire exceeds 
l-EmaxP/l-Ebl 2 > 2000, where E is the amplitude of the 
incident wave, as can be seen in Fig. [5] We note that 
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FIG. 5: Maximum field enhancement, |£7|/|i5o|, in the 
PIMNIM vs. wavelength. It is noted that the enhancement 
is maximized at the absorption peak, and has no structure 
corresponding to the electric resonance at A = 1380 nm. 



calculating the maximum field in the simulation domain 
can be imprecise, because mesh irregularities can lead to 
mesh-dependent spikes at the corners and edges of the 
structure. In order to calculate the field enhancement 
while avoiding such numerical artifacts, the results pre- 
sented in Fig. [5] were calculated by taking the maximum 
of \E\ interpolated on a rectangular grid of 1 x 1 x 1 nm 
and then smoothed by nearest-neighbor smoothing. This 
method was found to be robust under variations in mesh 
structure and density. Interestingly, the field enhance- 
ment is maximized at the absorption peak, at A = 1.5 /zm, 
yet shows no structure near the electric resonance, at 
A = 1.38 /im. This can be explained by the very strong 
reflection at the electric resonance, which prevents the 
incident fields from penetrating into structure. 

In addition to the enhancement near the corners and 
edges of the structure, fairly large field enhancement 
exists in a considerable volume between the cut wires. 
Fig. [5] presents \E\ at two y — z cross-sections, for dif- 
ferent values of x. In Fig. [6^a) the magnitude of the 
electric field \E\ is plotted at x = (in the middle of the 
PIMNIM). Intensity enhancement of \E\ 2 /\E \ 2 w 200 
is apparent in about two thirds of the region between 
the front and the back cut-wires. Fig. [UJb) presents the 
field closer the front cut-wire, at x = —8 nm. An av- 
erage intensity enhancement of |£?| 2 /|£'o| 2 s=s 400 is ap- 
parent in roughly the same area. If photon-counting de- 
tectors are integrated into the PIMNIM structure, then 
such intensity enhancement translates into proportion- 
ally enhanced absorption efficiency. The consequence of 
the large field enhancement is the desirable ultra-thin 
dimension of the absorber, in contrast to earlier calcula- 
tions 28 that demonstrated that thicker (about one wave- 
length) perfect absorbers can be developed using plas- 
monic spheres. 
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VI. CONCLUSIONS 



FIG. 6: (Color online) Field enhancement, |i?|/|_Bo| in the 
PIMNIM. Field is plotted on y - z slices at (a) x = and (b) 
x = — 8nm. Although the largest enhancement is near the 
edges and corners, a relatively large area experiences signifi- 
cant enhancement. 



In conclusion, we have demonstrated that an 
impedance-matched negative index metamaterial can act 
as a wavelength-selective wide-angle absorber of infra-red 
radiation. A specific implementation of such frequency- 
tunable PIMNIM based on plasmonic metamaterials is 
presented. Applications of the PIMNIM include in- 
frared imaging and coherent thermal sources. This work 
was supported by the ARO MURI W911NF-04-01-0203, 
AFOSR MURI FA9550-06-1-0279, and the NSF NIRT 
grant No. 0709323. 
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